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Abstract 

The realignment entanglement criterion of bipartite non- 
Gaussian states is obtained by applying the technique of 
functional analysis. The realignment criterion is given as 
one inequality in contrast to the infinitive number of in- 
equalities based on the moments. We give the necessary 
and sufficient condition of inseparability for non-Gaussian 
states prepared by photon subtraction or addition from 
symmetric Gaussian states. The entanglement criterion of 
non-Gaussian states evolved in thermal noise and ampli- 
tude damping environment is also obtained. 
PACS number(s): 03.67.Mn; 03.65.Ud 

Introduction — The inseparabilty of a non-Gaussian state 
is a longer standing problem in continuous variable (CV) 
quantum information theory. For the entanglement of 
a quantum state, various criteria had been developed, 
among them are Peres-Horodecki positive partial trans- 
pose (PPT) criterion [T], the uncertainty criterion [2] [3J, 
the entropy criterion [?] and the cross norm and realign- 
ment criterion [5] [BJ. The entanglement criteria initial 
developed for the case of discrete variable systems are also 
applicable to continuous variable systems. The entangle- 
ment criteria had been obtained for Gaussian states [7] [8] 
[5] [TU] which are completely characterized by the first and 
second moments. For 1 x l(each partite with one mode) 
[7] [S]and 1 x N llj Gaussian states, the entanglement 
condition is necessary and sufficient. It ceases to be suffi- 
cient for all higher mode systems [11]. The inseparability 
of non-Gaussian states are [T3] [12] Q3] [IS] [S] under in- 
tensify researches after the success for Gaussian states. A 
hierarchy of entanglement conditions were derived using 
high order moments [T2] . Although strict in the sense of 
PPT, they are not convenient for use except limited to low 
orders of moments. For each particular CV state, a large 
number of moments should be calculated when the orders 
of moments get higher and higher. The same problem ex- 
ists for Fock space PPT criterion [IS] . The newly founded 
entanglement criterion based on teleportation fidelity re- 
lies on the choice of test function [16] . Until now, to our 



best knowledge, all the entanglement conditions obtained 
for non- Gaussian states are sufficient but not necessary . 

We will use realignment criterion to study the entangle- 
ment of non-Gaussian states. The method is to write a 
non-Gaussian state as the functional derivative of a Gaus- 
sian state (Gaussian kernel). We first treat the realign- 
ment operation to the Gaussian kernel in Fock space, the 
operation then is expressed as a transform of the second 
moments. With the help of the generating functional and 
its characteristic function, the trace norm of the realigned 
non-Gaussian density matrix then is simplified. If the trace 
norm is greater than 1, the state is entangled. For some 
non-Gaussian states with experiment interests, the realign- 
ment entanglement criterion is necessary and sufficient. 
A non-Gaussian pure state will evolve to a non-Gaussian 
mixed state when it interacts with environment. There 
are debates [17] on the conclusions based on numerical 
calculations of the critical separable time of the evolved 
non-Gaussian states. We derive the strict formula of the 
realignment criterion to clarify the problem. 

Realignment entanglement criterion of Gaussian 
States — In Fock basis, a two mode Gaussian state with 
nullified first moments is written as 15 
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where s = (sx,S2),s' = (s'i,s' 2 ); o~i are Pauli matrices, 
and li is the 2x2 identity matrix. Here 7' = 7 + o\ ® 
Y , where 7 is the complex covariance matrix (CCM) of 
the state p defined by 7^ = ±Tr[p(ACiACj + ACjACi)], 
with C = (a\, a\, ax, a 2 ), AC = C-Tr(pC) and dj, ajare the 
annihilation and creation operators of the two modes. The 
matrix /? is 



(3= {(T3®l2W~ 1 (<T3®h)- 



(1) 



The first moments of the state can be removed by local 
operations, so Tr(pC) — is assumed. 
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The realignment operation transforms the state p to p R 
(p R is usually not a quantum state) such that [18] 



realignment entanglement criterion Tr y p R p R ^ — Trp R > 
1 is equivalent to (6 + ci) 2 — (b + c\) — c 2 < 0. That is 



Pktk-2 
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(2) 



The realignment operation leads to a chain of interchanges 
or substitutions. It can be expressed as the interchange of 
the subscripts: fc 2 wii, then can be expressed as the 
interchange of the parameters: s 2 s\, the operation can 
be further expressed with a change of the matrix j3 to /3r, 
we have 

P R = Z/3Z + I 2 ® o-i - CTi ® I 2 (3) 

where Z is a 4 x 4 matrix with entries Zn = Z23 = Z32 = 
Z44 = 1 and all the other entries are zeros. Denote 



7 - 1 = zy- 1 z' + / 2 



1 £71 + £71 



(4) 



(60 + ci + c 2 )(6 + ci - c 2 ) < 



(9) 



The second case is that 7r is not a valid CCM, then we try 
the operator p^ — p R H which may have all nonnegative 
singular values such that Tr\J p R p R ^ = Trp R , where 

n = J2 (-l) ni+n2 Imna) (nm 2 \ . 

n 1: n 2 



Then p R = J ^j^ p'n, wher e the CCM of the quantum 

state p' n is j Rn with 7^ = Z'^'^Z" + I 2 ® o\ + a\ ® 
7 2 ,with Z" = (c7 3 <g> I 2 )Z'. Then 



with Z' = (0-3 ® I 2 )Z(er3 <8> h)- Using local operations, 
the covariance matrix of a two mode Gaussian state can 
be brought to its standard form [7]- We assume 7 to be a 
standard CCM. A symmetric two mode Gaussian state is 
defined by Tr(pa\a\) — Tr{pa\a,2). Then for a symmetric 
Gaussian state 



ci 6 c 2 

ci c 2 6 

bo c 2 ci 

c 2 6 ci 



(5) 



where the real parameters C\ 
Tr(pa[a2) = Tr{po^ ' ■ 



- Tr{pa\a\) = Tr{paia 2 ), 
C2 = l r(pa[a2) = lr{pa 2 ai), b = Tr(pa\a\) + \. To 
simplify notations, we denote matrix 7 by its first row as 
7 = {0, ci, 60, C2}- Denote b = &o+~,then 7' = {0, c\, b, c 2 }. 
Hence 7 /_1 = Ai/ 2 <x>/ 2 + A 2 / 2 <8)C7i + A3(7i<g)/ 2 -|-A4C7i<g>C7i 
can be denoted as 



7 

with K x = 



= {^,^2,^3,^4}, 



(6) 
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{Kx,l- A 3 ,l + A 2 ,-A 4 }. 



The effect of appending II corresponds to reversing both 
signs of ciand c 2 . The realignment entanglement criterion 
is 

(&0-C1-C2X60-C1 + C2) < -.. (10) 



The condition for the entanglement of a symmetric Gaus- 
sian state is that either © or (fl"U|) should be fulfilled. The 
result of realignment entanglement criterion is exactly the 
same obtained previously for a symmetric Gaussian state 
by other criteria [7] [8] [T2]. 

The original realignment criterion for entanglement is 
Tr^J p R ^ p R > 1. In general case, we may not have 



Tr^J p R ^p R 



Trp R . Note that p R has the singular 



value decomposition p R = UAV\ with unitary opera- 
tors U, V and diagonal operator A = diag{Ai, A 2 , . . .}. 
So Try/p R ip R = Tr\A\ . Let U' = V^U be another uni- 
tary operator, we have Trp R = Tr(U'A) = £\ (7^-Aj < 
\ Ujj \ \Aj\ < \Aj\, the last inequality is due to the 
fact that \Ujj\ < 1 for any unitary operator U' . Hence if 



have Z'i- l Z' = {Kx, -A 3 , -A 2 , A 4 }, hence 



in 1 



{Ai,l- A 3 ,l- A 2 ,A 4 }. 



(J) 



2). 54 = fr(^ - c|- 6 2 )s where A 2 = det( 7 '). We Tr()R ig haye TrpR < Tr y^^ So if Trp* > 1 

then TrvW > Trp R > 1. 

For a non-symmetric state with standard covariance 
matrix, the inseparable criterion can also be written in 
the form of ([9]) or (| 10[) with 6q replaced by ^/bxb 2 , and 
6i = Tr{pa\a,i) + § fori = 1,2. 

Realignment inseparability of N on- Gaussian State — For 
the inseparability of non-Gaussian states, we introduce the 
method of generating functional. Let 



Notice that p R may not be a quantum state, however it 
is proportional to a quantum Gaussian state p' . There are 
two cases. The first case is that 7^ is a valid CCM, then 
p' is a quantum Gaussian state with 7^ being its CCM, 

and p R - 



d^TTF P ' so we S et 



Trp R 



/ det( 7 ^) _ / det(Y-i) 
det( 7 ') V/det^ 1 )- 



W = e eat e« a pe" at e Ca , 



(11) 



(8) 



A direct calculation shows that det(7^ 1 ) = det(7' 1; 



Q, where fl = 4[(A 2 
Further we have $1 = 



A3) 2 - (ATi 



A- 4 ) 2 ](l 



A 2 - A 3 ). 



^[(6 + ci) 2 -{b + cx)-cl\. Thus the 



where e, ^, 77, are two dimensional real parameters, we 
denote £idi + ^ 2 a 2 as fa and so on, p is the initial two 
mode Gaussian state with CCM 7. Then non-Gaussian 
state pn can be obtained from functional W by deriva- 
tions, pm — OW, where O is the operation of derivations 
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on the parameters e, £, 77, £ and later setting all these pa- 
rameters to be zeros. In Fock basis, 



where r 



4[(6o+ci) 



We have 



W klk2 



z v6 T +h(S,vMC,v) T d k, d k 2 g mi QP R = -f[(b - a T + ((&o + ci)r T + f(r + l) 2 



/mi 



v //ci!fc 2 !mi!m 2 !det( 7 ') dsj 1 0a£ 2 5s'/ 



<'M> 2( s ' s ')[^i ® ^2 +/3](s,s') T 



• exp { ( s , s ') [/?(e,^) T + (ry + e, C + C) T } L. 

The realignment is the interchange of k 2 4=> mi such 
that the realigned functional operator W has entries 
w k 1 k 2 , mim2 = W kimiMm2 . The operator W R can be put 
in the form of 



the new parameters are (£ R ,r] R ) = (£,r])Z; (e R ,( R ) = 
(e,()Z. Denote the realignment of the non-Gaussian state 
as p R = OW R . The realignment entanglement criterion 
is Trp R — OTrW R > 1. The characteristic function of 
the operator W is x W {^) = Tr[W V(p)], where V(p) = 
exp(/ia^ — /x*a) is the displacement operator, with p being 
a two dimensional complex parameter. After a lengthy 
calculation, we obtain 



for photon subtraction and addition. If the Gaussian ker- 
nel is on the separable boundary, namely, r = 1, then 
OP R = 1. This indicates that p T are separable. In fact, 
=jf=the Gaussian kernel p is separable thus it is proba- 
bility mixture of product states, the definitions of p± 
show that p± are also probability mixtures of product 
states. We confirm that if the Gaussian kernel is separa- 
ble, /9=p are separable too. If the Gaussian kernel state p 

is entangled, namely, yj ^ (y ) " = V^" > 1, we can show 
OP R > 1. Notice thatOP B | T=1 = 1, we need to show 
OP R \ ^ > OP R \ , , this is true from (115]). The neces- 

It>1 It— 1 7 ' ' 

sary and sufficient entanglement condition of p^\s exactly 
the same as that of its Gaussian Kernel state. Consider 
states prepared by subtracting or adding m photons from 
each mode of the standard symmetric Gaussian state p. 



The states prepared are p- m — c_ m a™a™ pa[ 



fm fm 



ax 



and 



P+m = c _ m c^™ a^ 2 m pa™ a™ . The states are separable when 
their Gaussian kernels are separable by the same reasoning 
for p± . Assume a function 



Pexp[--(p, p*)j(p, p*) T ] exp{(-/x, p*) 



(12) 



f(m,V) = 



d 1 



where 



c=0 



r 1 7 
exp[ — —KVK 



P = cxp [--(e+r ] ^+Or 1 (e+V^+0 T -Cv T -ee-ve}- 

Similarly, we can obtain the characteristic function 
X WR {IJ) of the operator W R , then TrW R = X WR (0) = 

£$^. with 

i = exphi^ + ^^ + CWC^ + ^^ + C^ 



for a matrix V — {Vi, V 2 , V3, V4}, where K — (fti, . . . , K4). 
Notice that Trp N = TrOW = Ox(0) = OP = 1, we have 



op r = /KTflT s^i® /a) 



/(m,7=F 9 o-i O J 2 ) 



(16) 



R c RT-\ 



The realignment inseparable criterion for a non-Gaussian 
state is 



/ det(7k) 0pfl > L 



(13) 



V det( 7 ') 

The first example is the non-Gaussian state prepared 
by photon subtracting from or adding to a symmetric 
Gaussian state. For a photon subtracted state p- = 
C-a\a 2 pa\a\ or a photon added state p+ = c + a\a 2 paia 2 „ 
where c T = [(&oT|) 2 +Ci+c|] _1 are the normalizations. A 
simple technique to calculate j' R is Yr = H^Ir 1 H) ~ X H, 
where H — I2®H2, with H 2 being 2x2 Hadamard matrix. 
Thus 

i'r = \{c2{t- 1),(6 + ci)t- (60 -ci), 

(60 + ci)r + (6 - ci) + 1, c 2 (r + 1)}, (14) 



for states pzp m . For the states p^ 2 , we have 

P f(2,V) = [V? + 2{V 2 2 + + V%)) 2 + Z2V 1 V 2 VzV i 
+8{V 2 2 V? + V 2 2 Vl + V^Vl). 

Numerical results show that OP R > 1 for t > 1. Hence 
the necessary and sufficient condition of inseparability for 
states p T2 is ^ too. 

The second example is a non-Gaussian state pn = PPi + 
(1 — p)p 2 , the probability mixture of an entangled Gaussian 
state p\ with a separable Gaussian state p 2 . For simplicity, 
consider the case of symmetric Gaussian state with c 2 = 0, 
namely the two mode squeezed thermal state (TMST). 
Denote w — bo + c\ — |, then the entanglement criterion 
is simply w < 0. Thus W\ 6 (— 1,0) for p\ and w 2 > 
for p 2 . The inseparable criteria of the second moments [5] 
, the Fock space [15] and the realignment for p^ are 

pwi + (1 - p)w 2 < 0, 
pwi + (1 - p)ui 2 + wiw 2 < 0, 
pwi + (1 — p)w 2 + 2w\w 2 < 0, 
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respectively. Hence realignment criterion is the best 
among the three. 

Inseparability of time evolution states — For the evolved 
state, we will use subscript t in the following to specify 
its operator and parameters, while the initial state is not 
with such a subscript. In an amplitude damping and ther- 
mal noise channel, the time evolution of the characteristic 
function for any quantum state is [19] 



Xt(p) = X(jue" 



: )exp[-(n+-)(l 



-rt\ 



l/'l 



(17) 



We have used the fact that p R is not a quantum state, 
however it is proportional to a quantum state p t with CCM 
jRt- Here 7^ is the realignment transformation of the time 
dependent CCM j t - The characteristic function oiW R can 
be obtained accordingly following the process from Eq. pTj) 
to Eq. ([T2|) . however, as we will show below, what we need 



(21) 



to know is Tr(W t R ) = xf R (0), which is 



xV R (o) 



PPT 



' det (7m) p fl 



det( 7t ') 



where T is the amplitude damping coefficient and n is the 
thermal noise, we assume the two mode undergo the same 
environment. Notice that X W (fi) m ([12} is proportional to 
the characteristic function of a two mode Gaussian state 
with first moments. The time evolution of X (t 1 ) then is 

xY(fi) = PcM-l(^^*ht{^,^*) T ]exp{(-fi,n*) 



where 



DH _ 



exp[-i( e f ■ + 



.RlcR\T 



■W'iet+ritAt+Ctf + ttuVtf}} (18) 



Where 



7t 
f% 



-rt 



-rt 



)<ti ® h 



= {<rz®h){~ft + vi®-£-) X (o- 3 ®I 2 ). 



The equation e 

Vt^t+(t) T 
tion 



Denote 
P = 



(?t; 7 7t) T leads to the parameter transforma- 



with p Rt = (a 3 ® hWm)' 1 (0-3 0/2). 

The characteristic function of the initial state is 



X n(p) = Tr[OWV(p)] = OTr[WV(fx)} = X W >)■ 

The second equality is due to the fact that the trace opera- 
tion is on the creation and annihilation operators crand a, 
and O is about the parameters e, £, n, £. The evolution of 
the characteristic function is a transformation of p, which 
is independent of the operations on parameters e, £, 77, C, 
hence 

We further have the evolved non-Gaussian state 
~d 2 /j,i~ 



rt 



exp[--( £i + + Ct)A '(et + ?7 t ,6 + Ct) 5 



n 



Xm(fi)1> (-fi) 

x7(fi)v(-n) 



then the time evolved W operator is 



ow t , 



and the realignment of the non-Gaussian state p R t = 
OW R . If all of the singular values of p R t are non-negative, 



then Tr^p R m p R \ = Trp R t = OTrW t R = 0*^(0). The 
realignment entanglement criterion for a bipartite non- 
Gaussian state is 



Wi 



PP^e eta 



ita p t e nta e CtQ 



(19) 



where p t is a Gaussian state with CCM 7$. The realign- 
ment operation on time dependent functional Wt will lead 
to 



/det( 7 ' 



ni 1 



OPPr x P R > 1 



(22) 



where 



R 



= PP t 



= pp; 



V det( 7t ') 

If some of the singular values of p R t are negative, then 



/det( 7 ^ 



det( 7t ') 



jM e e^ e t?a p ' te «*a\ 



(20) 



TryPmPm > Tr PNt- We sti11 have (EH) as the sufficient 
condition of entanglement. It can be further written as 



R 



= It 



h 



; det( 7 ^) 



Oexp[ 



1 



fe K ,^) 



(it,t)t)Z\ 

(et,Ct)Z. 



V det( 7t ') " ~- r " 2 
where v =(e, ~C, »7> -£)> 

M' = 7 " + e- rt 7 "( 



M 
M' 



M' 
M' 



™]>1. 



(23) 



/— 1 ^7/ / ry/ f—l\ I 

it z iRt z it h) 



"-1 

7i + 7t 



t/-l 7.;/ 7„,"-l\„," 

^7ra^7t )7 i 
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with 

7" = 7' -0-10/2) 

it = it-°i®h, 

ikt = i Rt -oi®h- 

Inseparability of photon subtracted and added two mode 
squeezed vacuum states — The initial state of photon sub- 
tracted state is p$ — c s aia2pa\a\, where p = \ip) (ip\ , 



71=0 



is a two mode squeezed vacuum state with A = tanh(r), 
where r is the squeeze parameter. The normalization con- 
The state is anti-normally ordered 



stant is c s 
with operator 



(l-A 2 )* 

X 2 (l+X 2 ) 



O = 



d 4 



=ri=0 



The initial state of photon added state is pA 



c a a\a\paia2, with c a 
dered with 



O 



1+A 2 



d 4 



The state is normally or- 



° deide 2 dCidQ 2 



e=C=0 



The CCM of the two mode squeezed vacuum state is 7 = 
{°'-T^ ,T^-iO}, so ^ = {0,-^ 7 _J_^ ;0 }and 
i t = {0, -A, N, 0} where 



A 
N 



X 



1 - A 2 
A 2 



We have 



im 



= {0,1 



-Tt 



,-rt 



N 



,-rt 



A 



7V2 _ a 2 ' 1 TV 2 - A 2 



,0}, 



so 



/ det( 7 k t ) 



V det( 7 2(N-A)-r 
The detail calculation shows 



(24) 



1—1 1 1—1 ryl I 17I 1—1 

-It + 7t Z lm Z 7t 



II— 1 1 //— 1 ry II ry II— 1 
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1 



2A^ - 2A - 1 



{0,1,1,0} 



for the evolution of two mode squeezed vacuum state. For 
the photon added state, we have r\ = £ = 0, the matrix 



M 



1 



1 + A 



{o,- 



A 



1 



1-A 



1-A 



O.o}, 



with Q = ( 1+ a)[2(]v-A)-i] e r *- The inseparable criterion 
is reduced to 



2 (w Vmi+V ( t A a + « 2 + (rb - > 1 

(25) 

For the photon subtract state we have e = £ = 0, the 
matrix M' = ^{O,-^ - XQ, A(j^ - Q),0}. The 
inseparable criterion is 

(26) 

The entanglement conditions based on second moments [5] 
[T21 are 



Ar- A + 



(l-A)V rt 



1-A 4 

for photon added state and 



e- 1 1 < 1 



N -A 



A(l - A) 3 _ r( 
1-A 4 



e- 1 r < 1 



(27) 



(28) 



for photon subtracted state. For photon subtraction, the 
critical time separability is comparable for the two crite- 
ria. Realignment criterion is better than the second mo- 
ment criterion at the higher initial squeezing parameter 
side (A > 0.3). For photon addition, realignment criterion 
is better for almost all the parameters. 

Summary — We derive a simple formula of the realign- 
ment criterion for bipartite continuous variable bipartite 
system. The main contribution of the paper is to pro- 
pose generating functional and its characteristic function 
in calculating the trace of the realigned density operator 
for arbitrary two-mode non-Gaussian state such that the 
realignment criterion is applicable to continuous variable 
system. We reduce the trace to the product of two fac- 
tors. One factor is determined by the second moments of 
the Gaussian kernel of the non-Gaussian state. The other 
is a derivative of an exponential function. For symmetric 
Gaussian states, the realignment criterion gives the same 
condition as the other criteria. For photon added states 
and subtracted states prepared from symmetric Gaussian 
states, the separable criterion is necessary and sufficient. 
It is just the separable criterion of the source symmetric 
Gaussian states. The condition can be extended to bi- 
photon (possibly to multiple photon) subtraction or addi- 
tion case. The realignment criterion is better than second 
moment criterion and one photon Fock space criterion for 
detecting the entanglement of the mixture of entangled 
TMST and separable TMST. The complete formula of the 
trace is derived for any two mode continuous variable sys- 
tem evolved in thermal noise and amplitude damping en- 
vironment. 

XYC and LZ J thank the support of the National Natural 
Science Foundation of China (Grant No. 60972071). 



■5 



References 

[1] A. Peres, Phys. Rev. Lett. 77, 1413 (1996); M. 
Horodecki, P. Horodecki, R. Horodccki, Phys. Lett. 
A, 223, 1(1996). 

[2] H. F. Hofmann and S. Takeuchi, Phys. Rev. A, 68, 
032103 (2003). 

[3] O. Giihnc, Phys. Rev. Lett, 92, 117903 (2004). 

[4] S. Walborn et al. Phys. Rev. Lett. 103, 160505 (2009). 

[5] O. Rodolph, Lett. Math. Phys. 70, 57 (2004). 

[6] K. Chen, L. A. Wu, Quant. Inf. Comput., 3, 193 
(2003). 

[7] L. M. Duan, G. Giedke, J. I. Cirac, and P. Zoller, 
Phys. Rev. Lett., 84, 2722 (2000). 

[8] R. Simon, Phys. Rev. Lett.84, 2726 (2000) . 

[9] X-B Wang, M. Keiji, and T. Akihisa, Phys. Rev. Lett. 
,87, 137903(2001). 

[10] G. Giedke et al. Phys. Rev. Lett., 87, 167904 (2001). 

[11] R. F. Werner and M. M. Wolf , Phys. Rev. Lett. , 86, 
3658 (2001). 

[12] E. Shchukin and W. Vogel, Phys. Rev. Lett., 95, 
230502 (2005) . 

[13] M. Hillery and M. S. Zubairy, Phys. Rev. Lett., 96, 
050503 (2006). 

[14] G. S. Agarwal and A. Biswas, New J. Phys., 7, 211 
(2005). 

[15] X. Y. Chen, Phys. Rev. A, 76, 022309 (2007). 

[16] H. Nha, S-Y Lee, S-W Ji, and M. S. Kim, Phys. Rev. 
Lett.108, 030503 (2012). 

[17] M. Allcgra et al. Phys. Rev. Lett. 105, 100503 (2010); 
J. Lee, M. S. Kim and H. Nha, Phys. Rev. Lett.107, 

238901 (2011); M. Allegra et al. Phys. Rev. Lett.107, 

238902 (2011). 

[18] C. Lupo, P. Aniello and A. Scardicchio, J. Phys. A, 
41, 415301 (2008). 

[19] X. Y. Chen, J. Phys. B, 39, 4605 (2006). 



6 



